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ABSTRACT 



We present hydrodynamic simulations of a major merger of disk galaxies, 
and study the ISM dynamics and star formation properties. High spatial and 
mass resolutions of 12 pc and 4 x 10^ allow to resolve cold and turbulent 
gas clouds embedded in a warmer diffuse phase. We compare to lower resolution 
models, where the multiphase ISM is not resolved and is modeled as a relatively 
homogeneous and stable medium. While merger-driven bursts of star formation 
are generally attributed to large-scale gas inflows towards the nuclear regions, we 
show that once a realistic ISM is resolved, the dominant process is actually gas 
fragmentation into massive and dense clouds and rapid star formation therein. 
As a consequence, star formation is more efficient by a factor of up to ~ 10 and 
is also somewhat more extended, while the gas density probability distribution 
function (PDF) rapidly evolves towards very high densities. We thus propose 
that the actual mechanism of starburst triggering in galaxy collisions can only 
be captured at high spatial resolution and when the cooling of gas is modeled 
down to less than 10^ K. Not only does our model reproduce the properties of 
the Antennae system, but it also explains the "starburst mode" revealed recently 
in high-redshift mergers compared to quiescent disks. 

Subject headings: galaxies: formation — galaxies: interactions — galaxies: star- 
burst — stars: formation 



Introduction 



Galaxies formed a significant frac tion of their stars during violent interactions and merg- 



ers, as probed by o bservations (e.g.. lElbaz &: CesarskyI |2003| ) and numerical models (e.g 



Hopkins et al.ll2006l ). The s trongest bursts of star for r nation are produced by rnaior merg 
ers of equal-mass galaxies (IBarnes fc HernquistI Il99ll : iMihos fc HernquistI Il996l : ICox et al. 
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20061). while the star formation (SF) activity rapidly decreases with increa sing; mass ratio 
dCox et alj[2008[ IPi Matteo et al.lboosi iManthev et aPboOsi iKniermanI 120091 ) . Nevertheless, 
the overall contribution of mergers to the SF budget of gal axies remains uncertain and qui- 



escen t star formation in isolated systems could domina te ( jjogee et al.ll2009t iRobaina et al. 



20091 : iForster Schreiber et al.l l2009l : iDaddi et al.l 1201 Obi ). In general, merger- induced SF is 



assumed to be centrally-concentrated , and to fuel the formation of bulges and compact 
spheroids rather than extended disks (INaab fc Burkertll2003l : iBournaud et al.l 120041 ). 



The mechanism usually invoked to explain merger-induced starbursts is that the inter- 
action with a companion induces an asymetry in the gas response, such as strong spiral arms 
and extended tidal tails, and the gas subsequently undergoes gravity torques. Inside the 
corotation radius, usually a few kpc, these torques are negative and the gas flows inwards: 
the density rapidly increases in th e central regions, and so do es the star formation rate. The 
mechanism was first detailed by iBarnes fc Hernquist (Il99in . a nd the resulting SFR were 
quantified in simulations by, e.g.. lMihos fc HernquistI (Il996[ ) and lCox et al.l ( l2006l ): the SFR 
in a merging pair can increase by a factor of a few tens compared to the same galaxies in 
isolation, but the increase i s generally much more modest - typically a factor of a few units 
compared to isolated disks ( jPi Matteo et al.ll2008l ). 



A fundamental and general observational fact about interacting galaxies is that their 
star formation proceeds in giant molecular clouds and star clusters, which can be 10-100 



times more massive than in normal spiral ga. 



of Super Star Clusters (SSCs; IWhitmore et al. 



axies, leading in particular to the formation 



20071). A theoretical explanatio n is proposed 



to be an increased gas turbulence (lElmegreen et al.l l2000l : IStruck et al.l l2005l ) , so that the 
Jeans mass, which sets the typical mass of gas clouds, becomes larger while the free-fall 
time of these gas clouds decreases. This mechanism of massive gas clouds formation could 
trigger the SF activity of interacting galaxies, independently of the traditional central inflow 
mechanism. 

However, most existing models of galaxy mergers, in particular those studying the SF 
activity, do not resolve clustered star formation in dense cold gas clouds: star formation 
is instead treated as a relatively smooth process (at least at scales of 100-1000 pc), taking 
place in a relatively homogeneous ISM supported by thermal pressure instead of a cloudy 
ISM supported by turbulent motions. Star cluster formation is then only indirectly modeled 
using sub-grid recipes (e.g., iLi et al.l 120041 ). This is caused by a limited spatial and mass 
resolution, and/or the absence of model for gas cooling below ~ 10^ K in many cases. Only 
a few models of galaxy mergers can directly resolve cold gas clouds and clustered SF (e.g., 
Wetzstein et al.ll2007l : iBournaud et al.ll2008l ) but the properties of merger- induced SF were 
not studied. 
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In this paper, we present adaptive mesh refinement (AMR) hydro dynamic simulations 
of a major galaxy merger. Our models have a maximal spatial resolution of 12 pc, and a 
refinement strategy ensuring that gas can cool down to a few 100 K. Gas fragmentation 
into dense clouds and star formation therein can thus be directly captured, at least down to 
masses of 10^ Mq. A pressure fioor avoids artificial fragmentation. A realistic multiphase 
ISM with den se clouds embedded in a warmer phase naturally arises in similar models of 
disk galaxies dXasker fc Brvanlbood : lAgertz et allboOQal : iKim et al.ll2009l ). 

Using this model for a merger of two galaxies with the interaction orbit of the Antennae 
galaxies, we study the interaction-induced SF properties and compare to lower-resolution 
simulations with a smoother, warmer ISM. We show that the process of gas fragmentation 
into massive clouds and rapid star formation therein dominates the merger-induced activity, 
while gas infiows become less efficient when a clumpy multiphase ISM is modelled. These 
results suggest that merger-driven SF does not follow the processes revealed in lower reso- 
lution simulations. Consequences include a potentially stronger starburst, but also a more 
extended distribution of gas and SF. With these new properties of merger-induced SF, our 
model can explain the properties of the Antennae galaxies. We also propose an i nterpretation 



for the different SF efficie ncies observed in quiescent disks and active mergers (jPaddi et al. 



2010al : iGenzel et al.ll2010f ). 



2. Model and parameters 
2.1. Gas physics and star formation model 



We use the AMR code RAMSES to evolve the dark matter and the stellar component 
using a Particl e Mesh solver, and the gas component using a second-order Godunov scheme 
( iTeyssierl 120021 ) . Throughout our study, we model star formation with a Schmidt law: the 
local star formation rate is = e^pgi^/tg, where = ^/^TtJ^2Gp^ is the free-fall time 
computed at the gas density pgas. The efficiency of star formation is controlled by the 
parameter e* ~ 1 %. Star formation occurs only in dense enough regions (molecular clouds), 
defined by the gas density nn being greater than some threshold value n,,. These two 
main parameters are usually calibrated using observations of nearby galaxies and the so- 
called Kennicutt-Schmidt (KS) law. The real efficiency is high with a high threshold, but 
models with a limited re solution have to use a low threshold combined with a low effici ency 
(jWada &: NormanI 120071 ). which globally reproduces the same KS law (lElmegreenl |2002| ) . 



The global SFR in our model (like most others) can be computed directly by inte- 
grating the Schmidt law over the gas density PDF above the threshold n*. The problem 




Fig. 1. — Top: RGB rendering of the system shortly after the first pericenter passage (blue: 
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therefore boils dow n to predict the PDF evolution during the merger process. As s hown 
by several authors (lElmegreenl l2002l : lElmegreen &: Scald |2004J : IWada fc NormanI 120071 ) , the 
multiphase structure of the ISM is built up from complicated processes, involving radiative 
losses (ultraviolet and infrared line cooling, molecular and dust cooling) as well as various 
heating mechanisms (cosmic rays and UV heating, supernovae and stellar feedback) and of 
course self-gravity. Surprisingly, numerical experiments have shown universal properties for 
the gas density PDF in isolated galaxies, with log-normal or power law distribution shapes. 
This was explained as a fundamental property of isothermal (or dissipative) self-gravitating 
turbulence. In this paper, our goal is to resolve this supersonic turbulence to compute 
self-consistently the density PDF and the resulting SFR in the course of a galaxy collision. 

We use for that purpose a simple thermodynamical model mimicking gas cooling due to a 
detailed balance between atomic and fi ne structure cooling and UV radiation heating from a 
standard cosmic radiation background (IHaardt fc Madaulll996l ). Assuming solar metallicity, 
we have computed the equilibrium temperature as a function of gas density to define our 
polytropic equation of State (EoS) T = Teq(nH). For densities between riu = 10~^ H/cc 
and riH = 0.3 H/cc, our model has T^q ~ 10^ K, while above nn = 0.3 H/cc, we have 
Teq ~ 10^ X (nH/0.3)-^/2 X. In this derivation, we neglected self- shielding of the radiation by 
the gas: gas cooling at very high density may have been underestimated, but we have also 
neglected the effect of local radiation sources such as OB stars as additional heating sources. 
Although our thermal model appears rather uncertain, it provides a reasonable route of gas 
dissipation, maintaining the gas temperature to a realistic average value at a given density. 
This EoS-based model produced a log-normal gas density PDF in isolated galaxies similar 
to the results of complete cooling /heat ing calculations (see se ction [3l) and a realistic density 
power spectrum of ISM substructures ( iBournaud et al.l 120101 ). 



Another import ant ingredient is the thermal support added at small scales to avoid arti- 
ficial fragmentation ( iTruelove et al.lll997l ). To ensure that the Jeans length is always sampled 
by at least 4 cells, we add an artificial pres sure defined as Pjeans = ^Q^x'^GPgas/l'^- This tech- 
nique, introduced in a different context by lMachacek et al.l (l200ll ). efficien tly prevents the for- 
mation of spuriously fragm enting clumps in galaxy formation simulations (IRobertson &: Kravtsov 
20081 : lAgertz et al.ll2009bl ). From our EoS, we can compute the typical density at which this 
pressure floor dominates, namely njeans — 6 x (Aa;/100 pc)~^/^ H/cc and the equivalent gas 
temperature as Tjeans — 2500 x (Aa;/100 pc)^/^ K. This density corresponds to our minimum 
thermal Jeans mass, which defines our mass resolution mres = iVfjeans- 

We use a quasi-Lagrangian refinement strategy: each cell for which the mass exceeds mres 
is subdivided into 8 children cell, down to the maximum level of refinement. In order to study 
the convergence properties of our system and identify qualitative changes, we performed a 
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low resolution model with Ax = 96 pc, mres = 10^ Mq and a high resolution one with 
As = 12 pc, rrij-es = 4 x 10"^ Mq. Note that in the low resolution run, the gas cannot 
cool significantly below 10^ K, while in the high resolution case, it can reach a minimum 
temperature of Tjeans — 500 K at the Jeans density. 

We have fixed the star formation efficiency parameter to = 0.01, adjusting the star 
formation density threshold to = 0.1 H/cc (resp. n,, = 8 H/cc) for the low (resp. high) 
resolution run. This ensures that both simulations initially have the same initial SFR of 
~ 1 M0/yr per galaxy in isolated disks, in agreement with the KS law of local spirals. 
Since in both cases, ra^, is significantly below n jeans, the star forming part of the PDF is well 
sampled. 



2.2. Merger parameters and initial conditions 



We model a pair of galaxies in a box of size 200 kpc, wi th isolated and ou tfiow boundary 
conditions. Each galaxy is embedded in a live halo with a iHernquistI (Il993[ ) density profile 
with masses Mh,i = Mh,2 = 1-8 x 10^^ Mq and scale lengths equal to the truncation radii 
'"h,! = ''"h,2 = 30.8 kpc. The pre-existing stars are described as two exponential disks of 
masses M*_i = 2 = 3.6 x 10^° Mq, scale length r^^i = rd,2 = 4.4 kpc and truncation radius 
'^max,! = 22 kpc, rniax,2 = 13.2 kpc respecti vely. The disks s cale height are = 0.2r^. A 
central bulge with B/D=l is added, with a IHernquistI (Il990l ) profile with a = 4.4 kpc. The 
gas distribution follows the stellar disk profile with a total gas fraction of 10% in each galaxy. 
The total number of dark matter particles was set to = 8 x 10^ and the initial number 
of stars to A^^om,* = 6 x 10^ f or both the low and high resolution simulations. We use the 
hyperbolic orbit proposed by iRenaud et al.l ( l2008l ) to reproduce the Antennae system. Our 
AMR grid has a coarse level £min = 7 (i.e. 128^ cells) and a maximum level of refinement of 
^max = 11 (resp. £max = 14) for the low (resp. high) resolution simulation. 



3. ISM dynamics and star formation in a galaxy merger 

Figure |2] shows the star formation history of our merger at low and high resolution. 
At first pericenter passage {t ~ 250 Myr), the SFR rapidly rises in both runs but with a 
dramatic difference in amplitude. Around t = 450 Myr (second pericenter passage) and 
t = 500 Myr (final coalescence), the SFR increases again, now more significantly for the low 
resolution run, and it steadily declines after the merger. 

Although the pre-merger isolated disks have similar SFRs, the global SFE during the 
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Fig. 2. — Star formation history in our simulations for the low (red) and high (blue) resolution 
runs. Times for first pericenter, second pericenter and final merger are indicated by vertical 
dotted lines. 
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Fig. 3. — Top panel: face-on view of the gas column density at time t = 350 Myr for disk 1 
for the low (left) and high (right) resolution runs. Bottom panel: corresponding fractional 
star formation rate as a function of radius. 
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Fig. 4. — Star-formation-weighted gas density PDF at three different epochs, before 
(t=160 Myr) and during the merger (t=345, 550 Myr), for the high (top panel) and low 
(bottom panel) resolution simulations. In each case, the star formation density threshold is 
indicated by the dotted vertical line. 
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Fig. 5. — Time evolution of the two merging galaxies in the KS diagram, comparing the low 
and high resolution runs. The solid line is a fit to the average K S law for quiescent galaxies, 
while the dotted line is the KS law for starburst galaxies (from lDaddi et al.ll2010al ). 
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mergers is totally different depending on the resolution, i.e. depending on whether we model 
a smooth and warm ISM or a cloudy multiphase ISM at high resolution. SFRs are discrepant 
by a factor of ~10 after the first pericenter and ~5 in the final merger stages. This is due 
to a completely different ISM dynamics building up two different density PDFs between the 
low and the high resolution runs (Fig. H]) and different spatial distribution of the gas (Fig. |3]). 

The low-resolution model, with a relatively smooth and warm ISM, follows the tradi- 
tional mechanism of merger- induced starbursts, where tidal torques excite a strong m = 2 
mode and drive gas inflows below the corotation radius. The typical m = 2 gas response 
and growing central concentration are seen on Figure |3l the enhanced star formation takes 
place mostly in the central kpc (Fig. |3]). 

The gas evolution largely differs in the high resolution run with a cloudy multiphase 
ISM. An m = 2 mode excited by tidal torquing is still visible, but the clumpy gas builds 
a much more modest central density peak than the low-resolution model: star formation is 
now spread over 2 kpc (Fig. |3]). 

The gas response to the interaction is actually dominated by fragmentation in many 
dense clouds scattered along the spiral arms, with a classical "beads on a string" morphology 
(see Fig. |3]), with typical clump masses of 10^~^ M©. These high clump masses result from 
the gas turbulent mot ions, which are increased by the interaction, as already studied in 
Bournaud et al.l ( 120081 ). The gas velocity dispersion in our isolated disk at high resolution 
is typically 10-15 km s^^ but it increases by a factor of 2-3 in most regions, and locally 
by a factor of 5 or more, during the interaction. The Jeans length increases in similar 
proportions, so the typical mass of gas clouds formed by gravitational instabilities at fixed 
average density is increased by factors of typically 10 to 100. Such increased gas turbulence 



is observed in interacting galaxies fe.g. . lElmegreen et a 



of supermassive gas clouds studied by lElmegreen et al. 



fll993l ). 



1995 ) and the associated formation 



Using the t heory of gravitat ional instability in a cylinder ( IChandrasekhar fc Fermilll953 



Qstrikerl 11964 : lElmegreen! Il979l ). we can compute the spiral arm stability criterion q = 
(T^/2G/i where /i is the linear mass density along the arm. We found /i ~ 4000 Mq/pc 
on average in both simulations. At low resolution, the minimum sound speed always lies 
above 5 km/s, so that g > 1 everywhere and the arm remains stable. At high resolution, 
the minimum sound speed can be lower (around 2.5 km/s) so that q can be as low as 0.25 
in quiescent regions, and the spiral arm can fragment into clumps in a free-fall time. As a 
consequence, the gas density PDF strongly evolves towards very high densities in the high- 
resolution model (see Fig. H]), and this evolution if achieved rapidly after first pericenter 
passage, and with little changes in the later stages: the star formation rate follows this 
evolution driven by ISM fragmentation in massive and dense clouds. 



- 12 - 



The low resolution model produces an artifically stable ISM and does not resolve this 
process. The gas response is dominated by the gradual inflow of gas towards the nuclear 
region with a timescale of a few 10^ yr (while dense gas clumps formed in a few 10^ yr). 
The associated density PDF increases slowly and continuously towards higher and higher 
densities as a result of the increasing gas concentration. This is why this model produces 
slower and more centrally-concentrated star formation than the high-resolution model. Note 
that our high resolution model do not resolve the clumpy ISM in the pre-merger isolated 
discs. This could affect our global SFR calibration in the initial conditions, but the initial 
gas density PDF would look very similar, except for the highest density tail. 



4. Implications 

4.1. The Antennae and other local mergers 

Although merger-induced star formation in our model is not primarily driven by an 
inflow of gas, and less concentrated than in earlier models, it remains relatively concentrated 
near the center of the merging systems (Fig. 3): there is still a tidally-induced inflow, and the 
central regions are denser and more prone to star-forming instabilities. This is consistent with 
star formation in ULIRGS being in general centrally concentrated (assuming ULIRGs are 
mergers). Nevertheless, our model also explains that the interaction- induced star formation 
can also be, for a part, radially extended (see Fig. [3]). This can explain why a number of 
interacting galaxies actually show extended star forna ation with SSCs forr ning at several kp c 



from their center, such as Arp 140 flCuUen et al.ll2006[ ) and the Antennae (IWang et al.ll2004j ). 



As for the Antennae, the orbit of which is matched by our simulation, there is a general 
consensus that we are witnessing the merger close to the second pericenter passage, when 



the two disks are still well separated ( iRenaud et al.ll2008l ). This corresponds to an epoch 



close to 450 Myr in our simulation (Fig. [T]). Our low resolution models reaches SFRs around 
10 Mq yr~^ just after the second pericenter passage and during only 50 Myr: these properties 
are in broad agreement with recent SPH simulations by Karl et al. (2010) and with the 
observed SFR (Zhang et al. 2001). Karl et al. (2010) then proposed that this high SFR 
and the extended SF in the Antennae results from the system being observed just at the 
particular instant of overlap between the two disks. Our high resolution model, however, 
shows that high SFR around 20 Mq yr~^ and relatively extended SF can be produced during 
a longer period (300 Myr) and does not require the system to be observed at a particular 
and brief instant. 



- 13 - 



4.2. General star-formation laws and the starburst regime 

Observations suggest a dual law for star formation, wtiere the integrated gas consump- 
tion timescale (HsFR/^gas) is relatively low for quiescent star-forming disks, and higher for 

nergers), as pointed out independently by 



Daddi et al. 


f2010a 


) and 


Genzel et al. 


(2010) 



tions, we retrieved integrated properties such as half-light radii, total gas mass and total 
SFR at several instants. The low-resolution model, where the starburst is driven only by 
the central gas inflow, does not show the observed change in 'LsFR/^gas'- the SFR increases 
during the merger, but only in proportion corresponding to the increase in the global gas 
density S^a^, and this model remains close to the standard relation for isolated disks (Fig. E]). 
The high-resolution model has its starburst driven mostly by increased gas turbulence and 
fragmentation. The gas density increases mostly on small scales in dense clumps throughout 
the system: this process does not affect the total effective size of the gas component, so the 
observable global density S^as has only a modest increase. At the same time, the starburst 
is even stronger than in the low-resolution model. This processes brings our model in agree- 
ment with the "starburst vs. quiescent KS law" pointed out by observations, throughout the 
duration of the merging process. Clustered star formation in high-resolution merger models 
can also affect the flnal structure of the resulting early- type galaxies (Bois et al. 2010). 

We thus propose that these recent observations unveiled a "starburst regime" where the 
efliciency of star formation on small scales and at high densities is unchanged, but exacer- 
bated gas turbulence and fragmentation into massive clouds result in faster gas consumption 
and higher integrated SF efliciency. The adopted SF law inside the clouds is not the key 
ingredient in the interpretation, the main eflect being the rapid evolution of the density 
PDF. Diflerent star formation models can indeed lead to diflerent quantitative predictions, 
but the main qualitat ive change comes from resolving the high density tail of the PDF 



(IGovernato et al.l 120101 ). Previous models of galaxy mergers did not resolve ISM turbulence 
and clouds, and could not unveil the physical processes driving this starburst mode. The 
actual process of star formation in starbursting mergers cannot be captured with sub-grid 
models on scales larger than 100 pc, and requires clustered star formation in a multiphase 
ISM to be directly resolved. 



Acknowledgments 

This work was granted access to the HPC resources of CINES and CCRT under the 
allocations 2009-SAP2191 and 2010-GEN2192 made by GENCI. 



- 14 - 
REFERENCES 

Agertz, O., Lake, G., Teyssier, R., Moore, B., Mayer, L., & Romeo, A. B. 2009a, MNRAS, 
392, 294 

Agertz, O., Teyssier, R., & Moore, B. 2009b, MNRAS, 397, L64 
Barnes, J. E., & Hernquist, L. E. 1991, ApJ, 370, L65 
Bois, M., et al. 2010, MNRAS, 915 

Bournaud, P., Combes, P., & Jog, C. J. 2004, A&A, 418, L27 
Bournaud, P., Due, P., & Emsellem, E. 2008, MNRAS, 389, L8 
Bournaud, P., 2010, lAU Symposium 270, Barcelona 

Chandrasekhar, S., & Permi, E. 1953, ApJ, 118, 116 

Cox, T. J., Jonsson, P., Primack, J. R., & Somerville, R. S. 2006, MNRAS, 373, 1013 

Cox, T. J., Jonsson, P., Somerville, R. S., Primack, J. R., & Dekel, A. 2008, MNRAS, 384, 
386 

CuUen, H., Alexander, P., & Clemens, M. 2006, MNRAS, 366, 49 
Daddi, E., et al. 2010a, ApJ, 714, L122 
— . 2010b, ApJ, 713, 686 

Di Matteo, P., Bournaud, P., Martig, M., Combes, P., Melchior, A., & Semelin, B. 2008, 
A&A, 492, 31 

Elbaz, D., & Cesarsky, C. J. 2003, Science, 300, 270 

Elmegreen, B. G. 1979, ApJ, 231, 372 

— . 2002, ApJ, 577, 206 

Elmegreen, B. G., Kaufman, M., & Thomasson, M. 1993, ApJ, 412, 90 
Elmegreen, B. G., & Scab, J. 2004, ARAA, 42, 211 
Elmegreen, B. G., et al. 2000, AJ, 120, 630 



- 15 - 

Elmegreen, D. M., Kaufman, M., Brinks, E., Elmegreen, B. G., & Sundin, M. 1995, ApJ, 
453, 100 

Forster Schreiber, N. M., et al. 2009, ApJ, 706, 1364 

Genzel, R., et al. 2010, arXiv: 1003.51801 submitted to MNRAS 

Governato, P., et al. 2010, Nature, 463, 203 

Haardt, P., & Madau, P. 1996, ApJ, 461, 20 

Hernquist, L. 1990, ApJ, 356, 359 

— . 1993, ApJS, 86, 389 

Kim, J., Wise, J. H., & Abel, T. 2009, ApJL, 694, L123 

Hopkins, P. P., Somerville, R. S., Hernquist, L., Cox, T. J., Robertson, B., & Li, Y. 2006, 
ApJ, 652, 864 

Jogee, S., Miller, S. H., Penner, K., Skelton, R. E., & Conselice, C. J., 2009, ApJ, 697, 1971 

Karl, S. J., Naab, T., Johansson, P. H., Kotarba, H., Boily, C., Renaud, P., & Theis, C. 
2010, arXiv:1003.685, submitted to ApJL 

Knierman, K. A. 2009, arXiv:0909.2053l 

Li, Y., Mac Low, M., & Klessen, R. S. 2004, ApJ, 614, L29 

Machacek, M. E., Bryan, G. L., & Abel, T. 2001, ApJ, 548, 509 

Manthey, E., Hiittemeister, S., Aalto, S., Horellou, C., & Bjerkeli, P. 2008, A&A, 490, 975 
Mihos, J. C., & Hernquist, L. 1996, ApJ, 464, 641 
Naab, T., & Burkert, A. 2003, ApJ, 597, 893 
Ostriker, J. 1964, ApJ, 140, 1056 

Renaud, P., Boily, C. M., Pleck, J.-J., Naab, T., & Theis, C. 2008, MNRAS, 391, L98 

Robaina, A. R., et al. 2009, ApJ, 704, 324 

Robertson, B. E., & Kravtsov, A. V. 2008, ApJ, 680, 1083 



- 16 - 



Struck, C, Kaufman, M., Brinks, E., Thomasson, M., Elmegreen, B. G., & Elmegreen, D. M. 
2005, MNRAS, 364, 69 

Tasker, E. J., & Bryan, G. L. 2006, ApJ, 641, 878 

Teyssier, R. 2002, Astronomy and Astrophysics, 385, 337 

Truelove, J. K., Klein, R. I., McKee, C. P., HoUiman, J. H., Howell, L. H., & Greenough, 
J. A. 1997, ApJL, 489, L179 

Wada, K., & Norman, C. A. 2007, ApJ, 660, 276 

Wang, Z., et al. 2004, ApJS, 154, 193 

Wetzstein, M., Naab, T., & Burkert, A. 2007, MNRAS, 375, 805 
Whitmore, B. C., Chandar, R., & Fall, S. M. 2007, ApJ, 133, 1067 
Zhang, Q., Fall, S. M., & Whitmore, B. C. 2001, ApJ, 561, 727 



This preprint was prepared with the A AS macros v5.2. 



